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Abstract 

We reconsider the study of persistent currents in a disordered one-dimensional 
ring threaded by a magnetic flux using he one-band tight-binding model 
for a ring of A'^-sites with random site energies. The secular equation for the 
eigenenergies expressed in terms of transfer matrices in the site representation 
is solved exactly to second order in a perturbation theory for weak disorder 
and fluxes differing from half-integer multiples of the flux quantum 0o = hc/e. 
From the equilibrium currents associated with the one-electron eigenstates we 
derive closed analytic expressions for the disorder averaged persistent current 
for even and odd numbers, Ne, of electrons in the ground state. Explicit 
discussion for the half-filled band case confirms that the persistent current is 
periodic with a period (/jq, as in the absence of disorder, and that its ampli- 
tude is generally suppressed by the effect of the disorder. In comparison to 
previous results, based on an approximate analysis of the secular equation, the 
current suppression by disorder is strongly enhanced by a new periodic factor 
proportional to 1/ sin^ for ^ (integer) </>o/2. 

PACS numbers: 71.50. -Ft, 72.15.Rn, 72.15.-v, 72.20.-i 



I. INTRODUCTION 

An Aharonov-Bohm flux threading a metallic or semiconducting ring leads 
to a persistent equilibrium current even if the ring is disordered, provided it is of 
mesoscopic size smaller than the mean free path for inelastic scattering. This is a 
direct consequence of quantum coherence which is not destroyed by elastic impurity 
scattering. The persisitent current is periodic in the magnetic flux with a period 
given by the elementary flux quantum 0o = hc/e (with h the Planck constant, c the 
speed of light and -e the electron charge) . The existence of persistent currents was 
predicted for one-dimensional rings by BYttiker, Imry and Landauer[l] and has later 
been observed experimentally in real three-dimensional rings by several groups[2-4]. 
Detailed theoretical investigations of the persistent current have appeared for one- 
dimensional systems as well as for multichannel two- and three-dimensional rings 
[5-28]. Several reviews of the subject have also been published [29-31]. 

Many of these studies were motivated by the experiment of Levy et a/. [2] 
whose result for the persistent current averaged over an ensemble of 10'' disordered 
copper rings is several orders of magnitude larger than the averaged current calcu- 
lated for models taking only the scattering by random impurities into account[14,16]. 
This led to attempts of including, in addition, the effect of the Coulomb interac- 
tion between the diffusing electrons[15,13,20,23,26], which improved the agreement 
between measured and calculated persistent currents considerably. However, in 
other recent numerical studies it was found that Coulomb interactions suppress 
the persistent current in one-dimensional disordered systems, rather than enhancing 
it[21, 24,25]. 

The present paper concerns the study of persistent currents in a one-dimensional 
disordered tight binding ring threaded by a magnetic flux, in the absence of electron- 
electron interaction. We recall that in previous discussions, starting with the key 
paper of BYttiker et al.\l], the effect of the flux is described via the flux-modified 
boundary condition (with x the position along the one-dimensional ring of circum- 
ference L) 

'il){x -\- L) — e'^'^'^%l){x) (1) 

for the ground state wavefunction of the system in the absence of fiux[32,33]. This 
description of the flux has been demonstrated for the Schrsdinger equation in con- 
tinuous space where the vector potential may be eliminated by choosing a gauge in 
which the external magnetic field vanishes in the interior of the ring [32,33]. 

We consider a tight-binding system composed of one-orbital atomic sites 
of spacing a on a ring threaded by a magnetic flux in the direction perpendicular to 
the ring. This system is described by the set of difference equations 

-e''^o^^^^-e''^'^oip^_^ + eniPn = E(pn,n = 2,3,...,N-l , (2) 
-e 'I'o (p2 - e 'I'oipN + ei(pi ^ E(pi , (3) 

-e^t'oif^-e 't-o(pN_i + eN^N = E(pN , (4) 
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which correspond to the familiar Anderson model for a disordered ring modified 
by the flux of the external vector potential [6]. Here v^n is the amplitude of an 
eigenstate wavefunction at site n, E and En are the corresponding eigenvalue and 
random site energies (potentials) in units of minus a constant nearest-neighbour 
hopping parameter. 

For clarity's sake we now briefly rederive the boundary conditions (1) for 
the discrete ring described by equations (2-4). To this end we perform a gauge 
transformation 



V'„ = e*'^^o>„, n = l,2,... ,Ar , (5) 
to new amplitudes ipn, which transforms the system (0-|D into 



-1pn+l - = {E - en)i^n , U = 2, 3, . . . , N - 1 , (6) 



.^2-e-*'"^^;v = (^-£i)V'i , (7) 



-e'^'' fo ijji - ipN-i = {E - eN)il>N , (8) 

where the effect of the flux is transferred entirely to the bond linking the sites 1 and 
N. These sites are bounding the continuous sequence of sites n = 2, 3, . . . , — 1 
described by the equations (|) from which the flux is eliminated. We now observe 
that the system (6-8) is formally equivalent to 



-iJn+l - = {E - En) i^n, n = 1,2, N , (9) 

with the conditions (where (10. b) is the analog of (1) taken at the origin) 



i>N+i = e^'^i'o^pi , (lO.a) 

-ipN = e '^^ipo , (10. b) 

which describe a flux- free system deflned by the N equations (9) , involving + 2 
amplitudes ipn, in which the flux enters via the boundary conditions (10. a, b). Indeed, 
these boundary conditions ensure that the equations (9) for n = 1, A^ coincide with 
(7-8). In addition, one sees that rotating all amplitudes in (9) once around the 
ring {ipm i^m+Nyfn = 1;2, ... ,A^) Icavcs these equations unchanged under the 
boundary condition 



ijjm+N = e^'^'/'Olprn ■ (10. c) 



This shows that the flux modified (or twisted) boundary condition (1) is also valid 
for the discrete ring. We further recall that the similarity between the boundary 
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condition (|T]) and the property of Bloch functions in a periodic potential[l] (with 
playing the role of the momentum and the ring size L corresponding to the 
size of the unit cell) implies a flux periodicity with period (j)o of the eigenstates and 
eigenenergies, and hence of the equilibrium persistent current carried by an energy 
level E„, 



dE„ 

/„ = -c^ . (11) 

Cheung et al. [5] have studied the persistent current in the system described 
by equation (9) and (10. c) in terms of a two steps process: in the first step the 
flux-free ring is being cut open and deformed int o a linear chain whose Schrsdinger 
equation (9) is transformed to a tight-binding Bloch wave basis [5,35] and is solved 
in terms of a transfer matrix for amplitudes of Bloch waves travelling to the right 
and to the left of the chain; in the second step the wave amplitudes at the two ends 
of the chain are further matched by imposing the flux modified boundary condition 
(10. c). This leads to an exact secular equation for the eigenenergies of the ring 
involving only the fiux and the total amplitude transmission coefficient, t, of the 
disordered chain of the first step at energy E. Cheung et al. combine their secular 
equation with (11) to calculate the persistent current, by making simple ad hoc 
assumptions about the energy dependence of t and approximating the effect of the 
disorder in the averaged transmission coefficient in terms of the Landauer formula 
[5]. In the case of weak disorder they find a constant, fiux-independent, reduction of 
the amplitude of the persistent current with respect to its value for a perfect ring. 

From a first-principles point of view the constant reduction of the current 
amplitude due to the disorder obtained by Cheung et al. [5] is surprising. After all, 
the exact explicit expressions of the transmission- and reflection coefficients, t and r, 
are complicated functions of the random site energies ej and of energy (through the 
dependence on the Bloch wavenumber, k = arccos(— i?/2)) [5,35], from which one 
may suspect the existence of a non trivial interplay of disorder and flux in persistent 
currents in the exact eigenstates. 

The purpose of this paper is precisely to discuss an exact flrst principles an- 
alytic calculation of the eigenenergies and persistent currents for the ring described 
by (2-4) in the framework of a second order perturbation theory for weak disorder. 
In particular, we wish to study quantitatively the effect of the expected interplay 
of disorder and magnetic flux on the magnitude of the persistent current. On the 
other hand, a way of assessing the relative importance of electron interaction effects 
in persistent currents [15,13,20,21,23-26] in the disordered regime, is by comparing 
numerical calculations with results obtained for non-interacting electrons. Exact an- 
alytical results for the effect of disorder on the persistent current for non interacting 
electrons, such as those discussed in this paper, are clearly valuable in this context. 

The perturbation theory for the eigenenergies of the ring for weak disorder 
is conveniently developed by starting from a secular equation expressed in terms of 
transfer matrices for the Schrsdinger equation (2-4) [or (6-8)] in the site representa- 
tion. The secular equation for the eigenenergies of the ring as well as their explicit 
determination to second order in the random site energies are discussed in Sect. II. 
In Sect. Ill we use the results of Sect. II for calculating the persistent current av- 
eraged over the disorder for a system of Ng spinless particles in the ground state 
of the tight-binding system. Our results will thus be relevant to average currents 
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obtained from measurements for a large number of weakly disordered mesoscopic 
rings corresponding to different realizations of the disorder[2]. The discussion of our 
results for the typical low filling {Ne <C N) and half-filled band (A^e — y) cases, 
their comparison with previous work and some concluding remarks are presented in 
Sect. IV. 
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11. TRANSFER MATRIX PERTURBATION 
THEORY 



Our aim is to calculate the electronic eigenenergies of the tight-binding ring 
described by (0-|D to successive orders in small random fluctuations Ei of the site 
energies about a zero mean value. This amounts to choosing the energy level of 
the individual atomic sites as the zero of energy. We flnd that the method based 
on tranfer matrices in the site-space representation[34,35] is ideally suited for this 
purpose since these matrices depend linearly on the energies Ei. 

Deflning the unimodular transfer matrices 

P-(-^^-'"^ a = ^^ n=12 N (12) 



we have, from 



^"+M=e— PjJ" ) ,n = 2,3,...,iV-l , (13) 



^^^-e-A !M , (14) 



^ - --P^ ( ) . (15) 



By combining the amplitude vector of components (p^ and ^n-i-, obtained by itera- 
ting ([131) in terms of the vector with components v?2 and with ([T^ITSl) , we readily 



obtain the following exact equation for the eigenvalues of the ring: 



N 



detm '^''^o JJ A - 1 ) = , (16) 
which simplifles to 



n=l 



N 



tr 



W PA =2cos27r^ , (17) 



since the determinant of the product of unimodular matrices, Y{n=i ^n, is equal to 
unity. 

For the purpose of our perturbation treatment for the eigenvalues for weak di- 
sorder we flrst discuss the solutions of (0) in the absence of disorder. 
In this case 



^A=p^withP=^ £ ) , (18) 

n=l ^ 
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where denotes eigenenergies of the perfect ring. In order to determine the ele- 
ments of the matrix we diagonalize P by means of a similarity transformation. 
The eigenvalues of P are 

Ai,2 = ^ [-E^ ± VE^^-A^ , (19) 
which may be rewritten in the form 

Ai,2 = e±^^ , (20) 

with 

^o = -2coss . (21) 

The similarity transformation which diagonalizes P is defined by the non-unitary 
eigenvector matrix 

/ i(s-a) -i(s+a)\ 6*° / 1 — 

Using the fact that {u^^PJJ^ = Xi6ij, i = 1,2, we then obtain from (|2lD and 

(IH), 



p"' = u(u-^PuYu-^ = J—(^^'^.^^^'^^-'^ -e *"sinms\ ^ 
V / sins \ e*" sin ms — sin(m — l)sy 

In the absence of disorder the eigenvalue equation (|T7|) then reduces to 



cosiVs = cos— — , (24) 

00 

which leads to the energy levels 



E°(0) = -2cos^(A: + ^) , (25) 



using ([2l|) , and to the corresponding equilibrium currents (equation ([Til) ) 



/° = -^sin^(fc + -p) , (26) 



Nh N 



^0 



where k = 0, ±1, ±2, .... These results coincide with those obtained earlier by 
Cheung et al.[b] for the perfect tight-binding ring. For levels with \k\ N the 
energies ([25|) as a function of 0/0o grow quadratically with 0/0o away from minima 
at 0/00 = —k, k = 0, ±1,.... For such levels therefore, the energy spectrum 
is similar to the familiar free particle spectrum composed of intersecting parabolas 
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whose ordinates increase as a function of </>/</>o away from their respective vertices 
at 0, ±1, ±2,... [5]. 

We proceed with the solution of (|T7|) for weak random fluctuations of the site 
energies where we expand the equation to successive orders in the Si and obtain the 
corresponding first and second orders corrections, E^^^ and E'j^\ for the perturbed 
energy levels, 

Ek = El + + ) + ■ ■ ■ , = 0, ±1, . . . . (27) 
The transfer matrices P„ are split into unperturbed and perturbed parts: 

P„ = P + V; = P-(AE-£„) °^ , (28) 



where 



AE = ^(1) + ^(2) + . . . , (29) 



and the matrix product, Y[n=i ^n, is expanded to quadratic order in the quantities 
AE-En. This yields 



N N 

SN-m 



n=l m=l 

N n-1 



71=2 711=1 



By inserting this expression in ([ITD , using (|23D, and equating the two sides of the 
equation order by order we obtain, after some straightforward algebra, 



N 

~ N 

71=1 



E'^'-^t^r, , (31) 



and 



N 71-1 



sm qk sm Nqk ^ ^ ' 



n=2 771=1 

sin(n — m)qk. sin(A^ — n + m)qk , (32) 



where 



5^ = 1^ ( ^ + i- ) • (33) 



'Jo 



Clearly, the weak disorder perturbation theory breaks down at flux values equal to 
half-interger multiples of the flux quantum 0o where ([3^) diverges. Such divergencies 
are familiar[19] in other perturbative studies of persistent currents[13,14,16]. 
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III. AVERAGED PERSISTENT CURRENT 



The random eigenenergy corrections (|3lD and (IS^) are easiest to analyse 
through low order moments. Here we shall consider only the averaged eigenenergies 
and the corresponding averaged persistent currents. These are relevant to experi- 
ments obtaining the persistent currents in a large number of identical mesoscopic 
rings corresponding to different realizations of the disorder [2]. In averaging the 
eigenenergies ( [27| ) we assume a site energy disorder in the form of gaussian white 
noise with a correlation 



From (|3l|) and (^) we thus obtain (E^^^) = and 



(34) 



N n-1 



A^^ sin Qk sin Nqj, 



sin(?T, — m)qk sin(A^ — n + m)qk ■ (35) 



n=2 m=l 

By performing the finite summations over sites in (R^) we finally get 



4 sin qk 



1 + ^ j cot Nqk - ^ cot qk 



k = 0, ±1, ±2,... . (36) 



From ([TT|), (|25D and (^) we then obtain the averaged current from the level 



(4) 



e I . Nel 

2 sm qk - — 

4 sm qk 




1 1 27r0 \ 1 / 2 N 

— ^ + ^cot--cot g.J +^(l + 2cot g,) 

</>0 



(37) 



which reduces to 



(4) 



Nh 



2 sin — ( A; + — 



4 sin 



2 27r^ 
<^0 



k = 0, ±1, 



, (38) 



by dropping terms of order and in the square brackett of (^7|). 

In order to obtain the total averaged persistent current (/) for a system of 
electrons in the ground state at a given 0/0o we must sum the contributions from 
the lowest occupied levels at the considered fiux value. For simplicity we ignore 
spin and assume that a level may be occupied only by one electron. As is well 
known the current is different for even and for odd numbers of electrons[5]. We first 
concentrate on the contribution to the total current from the zeroth order levels 
corresponding to the non-disordered system. For odd Ne and for values of flux in 
the interval ~| < ^ < | this contribution, is obtained by summing the one- 
electron currents (UB) from k = — '^"'^ to k = ^"'^ . The fact that these values of k 
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correspond to the Nf, lowest levels of the spectrum in the considered range of fluxes 
is easily verified e.g. in the free particle-like limit of ( pSj ) where the eigenernergies 
reduce to intersecting parabolas whose vertices on the 0-axis correspond to integer 
values of 0/0o- One finds in this case [5] 



(Ne~l)/2 2l!I± 



/(«)= 4" = -/o-^ -2-£/„ . (39) 

where 



Sm -r7 00 



NdEl 
h dk 



/o = ^, with t;^ = - sin ^ . (40) 

The meaning of the velocity Vp is clear: its value coincides with the quantity 
which is the band velocity at an energy just above the highest occupied 

k=Ne/2 

level, that is a velocity close to the velocity at the fermi level. It may also be simply 
related to a fermi velocity as follows. Consider a non-disordered linear tight-binding 
chain whose electronic states are described by the equation —V{ipn+i+^Pn-i) = E ipn 
for a system of sites. In the continuum limit 

(a 0) this equation becomes (with x = no) —a'^VdPip{x)/dx^ = {E — 2V)(p{x), 
which is the Schrsdinger equation for a free particle of mass m = h^/{2Va^). In this 
limit the fermi velocity for a system of Nf. electrons is 
Vp = hkp/n = ^^jf- (with kp = NeTT/{Na)) the fermi wavenumber). This 
shows that fi;' is the limiting value of the band velocity in (|40D for low band filling, 
NeTT <^ A^. This discussion shows that /q is the current associated with a single elec- 
tron of velocity Vp at the fermi level (or of fermi velocity vp for low band filling). This 
is the well-known manifestation of the strong compensation of the contributions to 
the persistent current from successive occupied levels of the system, with the result 
that the total current is approximately equal to the current from the highest occupied 
level. 

Similarly, from the structure of the level spectrum (p5|), it follows that for 
even states for 1 < 0/0o < the N^. lowest occupied states correspond to 
k = 0, ±1, ±2, . . . , ± (^ — l) , — In this case one thus obtains[5] 



sm 

/(°) = -/n — 



Af Uo 



-('^-1) ^0 . (41) 



Obviously the strong compensation of contributions from adjacent levels to 
the pure system current shown by ( ^0| - |4lD does not occur when summing the effect 
of the disorder in ( p^ ) for the successive occupied levels. Thus, unlike for the con- 
tribution from the non-disordered energy levels, we expect to obtain a reasonable 
estimate of the effect of the disorder on the averaged current by approximating the 
discrete sum over occupied levels by an integral. More precisely, we shall use the 
Euler-MacLaurin summation formula[36]: 
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n-1 
k=l 



f{k)dk - - [/(O) + f{n)] 



(42) 



where we have included only the leading integrated terms [36]. 

For odd N^. and ~| < ^ < | where the summation over occupied levels runs 
over = 0, ±1, . . . , ± (^^^Y^) we obtain, by applying (|42D, 



E 



smr-j 



smr_ 



(1 + cosr_,.)(l — cosr_ 
— In 

All (1 — cosr+)(l + cosr_ 



(43) 



where 



(44) 



In the case of even N^. , for < 



< 1, the summation over k is over the values 



A; = 0, ±1, . . . , ± - 1) and and we find 



1 1 
2 V sin s+ sin s_ 



(1 + cos s+)(l — cos s_ 
47r (1 — cos s+)(l + cos s_ 



(45) 



where 



vr 



s± 



(46) 



After some simplification of the disorder contributions the expressions for the total 
averaged current, (/) = J^ki^k) given by (|38D , (|3^) and (|4^ - |4^ for odd Ne and by 

Na I 



P^), ( |4l|) and ( ^5] - ^^ ) for even Ai'e, respectively, reduce finally to (with kp 



(I) = -h 



sm 



N<Po 



sm^ 



+ 



sin^ sin /ci?a 



cos(A;ira + -^) sin jy^ ^sin (fc^^a + — sin 



sm 



2 _27r. 



^1 



^ - sin^ (fciT'tt + ;|) 27r \^ sin (fci^a + ;|) + sin 



for odd iVp 




(47) 



and 
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IV. DISCUSSION AND CONCLUDING REMARKS 



In discussing the averaged persistent currents (|47|) and (|48|) we consider suc- 
cessively two limits for the so-called band filling factor N^/N: the case of low 
filling,^ < 1, A/^e > 1 , and the half-filled band case, = N/2, for even 
A^. For 71 Ne N , one easily verifies that the disorder corrections in (|T^-|48|) are of 



opposite signs to those of the perfect systems currents i.e. the disorder suppresses 
the persistent current. On the other hand, in the half-filled band case, = y, we 



obtain successively from and 



(/) = - 




Nel 



i sm 



2 2TT(t> 



and 



10 1 

for odd Ne, -77 < — < 

2 00 2 



7^0 , (49) 



{!) = - 




for even A^^e, < — < 1, — ^ - 

00 00 2 



(50) 



to leading order in The equation ( ^ - |50|) show that for the half-filled band 

cases the effect of the disorder is simply to multiply the perfect system current by 
a periodic fiux dependent reduction factor of period 0o/2. This reduction factor 
does not alter the periodicity with period 0o of the original perfect system current 
discussed in[5-7]. We recall that the reduction of the current due to the disorder has 
a simple interpretation in terms of level repulsion caused by the random impurity 
potential. Consider, for simplicity, the limit where the conducting electrons are free 
and the energy level spectrum as a function of 0/0o is composed of intersecting 
parabolas. As is well-known, the one-body impurity potential causes level repulsion 
which leads to the replacement of the free parabola crossings by avoided crossings. 
The latter are more fiat than the original true crossings, which leads to a reduction 
of the currents {Ik) = —cd{Ek)/d(j) near them, caused by the interaction with the 
impurities. 

The expression (^9| - ^ ) for weak randomness differ remarkably from the cor- 
responding results of by Cheung et al.{5]. While our results are obtained from an 
exact solution of the secular equation (17) to second order in the random site en- 
ergies for weak randomness, those of Cheung et al. [5] follow from an approximate 
study of their exact eigenvalue equation 



Re 



cos 2n- 



(51) 



obtained from (9) and (10. c) after transforming to a Bloch wave representation as 
discussed in Sect. I. We recall that t is the complex transmission coefficient at energy 
E across a disordered linear chain of length Na obtained by cutting the ring open. 
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From simple but crude assumptions about the energy dependence of t Cheung et al. 
derive an expression for the current from a level Ek which may be approximated by 
Ik ~ /o|t| [5]. For a half-filled band, assuming that the total current in the presence 
of disorder is approximately equal to — (20/0o — 1) times the current carried by the 
last occupied level they obtain [37] 



(/) = -f^-l) /o(|t|) , (52) 



^0 



where (|t|) is the averaged transmission coefficient measured at the fermi level which 
may be related to the localization leength by using the Landauer formula. For weak 
disorder (|t|) is thus found to decay as (in our notations [38]), with a constant 
coefficient, from the value one for the ordered system (See Sect. III.C of [5]). In 
contrast the exact weak disorder effect in (^), which is also proportional to Ne^, 



varies strongly with the flux, as 1/ sin^ The strong flux dependence of the 

disorder effect in the persistent current ([49| - |50D shows that the approximations used 
by Cheung et al. [5] in solving ([5l[) are insufficient, at least for the case of weak 
disorder. 

In discussing the magnitude of the effect of the disorder in (^9[-^) relative 
to the perfect system current, it is useful to recall the domain of validity of these 
perturbative expressions. This domain is determined by: 

1. the condition, \eo\ <C 1, of smallness of the amplitude fluctuations of the site 
energies relative to the hopping parameter in ([3fl) 



2. the inequality ^/ 2) ^ A, which expresses that the typical shift of energy 
levels due to disorder must be small compared to the spacing A of energy levels 
at zero flux at the fermi surface[5,7]. From (^) we have A = sinkpa, and 



by using (pT[) and ([3^ this condition reduces to 



ViV|eo|<47r , (53) 

in the half-filled band case. This implies that within the perturbative domain 
the disorder term in the brackett of ([50|) may be enhanced to values of order 
unity by the factor 1/sin^^, for fluxes sufficiently close to multiples of (j)o/2. 

For odd this leads to a strongly enhanced suppression of the current near = ±^ 
while for even A'^g a similar enhanced suppression occurs near = and near 
= ±00. Such large effects of a weak disorder may exist only for mesoscopic 
rings since e.g. for l^ol = 0.1 values of N compatible with (|53|) would be restricted 
to ~ 10^. 

We conclude by suggesting that important modifications of the amplitude of 
persistent currents in discrete rings due to the interplay of disorder and magnetic 
flux may also exist for more realistic multichannel disordered rings as well as for 
models incorporating Coulomb interaction. In view of the persisting disagreements 
between calculated and measured persistent currents in metallic rings (see Sect. EI) 
such effects would seem to deserve further detailed study. 
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